Superfluid Fermi gas in a ID optical lattice 
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We calculate the superfluid transition temperature for a two-component 3D Fermi gas in a ID 
tight optical lattice and discuss a dimensional crossover from the 3D to quasi-2D regime. For the 
geometry of finite size discs in the ID lattice, we find that even for a large number of atoms per 
disc, the critical effective tunneling rate for a quantum transition to the Mott insulator state can be 
large compared to the loss rate caused by three-body recombination. This allows the observation of 
the Mott transition, in contrast to the case of Bose-condensed gases in the same geometry. 
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The observation of a BCS superfluid transition remains 
a cfiallenging goal in the studies of ultracold Fermi gases. 
It was recently suggested that gases confined to low di- 
mensions are promising candidates for achieving super- 
fluidity as the confinement enhances interaction effects 
0. Adding a tunable periodic potential allows one to 
combine the benefit of the reduced dimensionality with 
the advantage to work with large yet coherent samples. 
In particular, it has recently been shown that the study of 
the center of mass oscillations of the cloud in a ID lattice 
plus a superimposed weak harmonicpotential allows one 
to probe the superfluid transition Importantly, the 
presence of periodic potential introduces a much richer 
physics related to a possibihty of observing a variety of 
quantum phase transitions 0]. 

In this Letter we obtain the BCS transition tempera- 
ture Tc for a two-component Fermi gas in a ID optical 
lattice, assuming that the Fermi energy is small com- 
pared to the interband gap and, hence, superfluid pair- 
ing occurs only in the lowest band. This requires us to 
reveal how the presence of the ID lattice renormalizes 
the effective coupling constant at the Fermi level. For 
the geometry of finite size discs in the ID lattice, we also 
discuss the possibility of achieving the Superfluid-Mott 
Insulator quantum transition by tuning the lattice depth 
above a critical value. In this peculiar phase, the gas is 
superfluid in each separate disc but the coherence along 
the lattice direction is completely lost. We show that for 
Fermi superfluids the critical effective tunneling rate can 
be large compared to the loss rate of all inelastic processes 
and therefore the Mott transition can be achieved. This 
result is a direct consequence of the Fermi statistics and 
is in marked contrast with the case of Bose-Einstein con- 
densates in the same geometry, where the Mott transition 
can be hardly observed as pointed out by Hadzibabic et 
al. 01, unless the number of atoms per disc is very small. 

We consider a two component atomic Fermi gas in the 
presence of a one dimensional (ID) optical potential 



Vopt = sEr sin qsz, 



(1) 



intensity of the laser beam. En — h^qg/2m is the recoil 
energy, with hqs being the Bragg momentum and m the 
atom mass. The potential has periodicity d = tt/qb 
along the z-axis. The weak attraction between atoms 
in different internal states is modeled by a s-wave pseu- 
dopotential U{r) — gS{r)dr{r-) with coupling constant 
g = ATrh?a/m, where a < is the 3D scattering length. 

We will discuss the situation where the laser intensity 
is sufficiently large (s > 5) and the Fermi energy e^^ is 
small compared to the interband gap eg . We thus confine 
ourselves to the lowest Bloch band where the physics is 
governed by the ratio of the Fermi energy to the band- 
width At, where t is the hopping rate between neighboring 
wells. For ep < 4i the Fermi surface is closed and the 
system retains a 3D behaviour, whereas in the case of 
> 4i the Fermi surface is open and the system under- 
goes a dimensional crossover. Hence, one has two distinct 
regimes: an anisotropic 3D regime {ep ^ t) and a quasi- 
2D regime (e^, 3> cf S> t). This is clearly different from 
the case of a 3D lattice 5] where the Fermi energy scales 
with the bandwidth and can therefore be much smaller 
than the corresponding value in free space for a given 
atom density. 

The mean field transition temperature is the high- 
est temperature at which the Gorkov equation for the 
gap parameter has a non-trivial solution This gives 
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(2^)3 Cqexp(Wr°) + l' 
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where ^e// is an effective coupling constant. The sym- 
bol P stands for Principal value and — ti^q^/2m + 
^i{Qz) — where q_L is the momentum in the direction 
perpendicular to the lattice, e(gz) is the band dispersion 
and fi ~ ep is the chemical potential. A straightforward 
integration of Eq. ^ yields 
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where s is a dimensionless parameter coming from the 



with 7 = 1.781 and = / (5($q)dq/(27r)3 being the 

density of states per internal state at the Fermi level. 
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The function F is defined as 



F ■ 



/!-^dg,ln(l-efe)/M)eO.-e(g,)) 



(4) 



where Q{x) is the unit-step function. 

The effective couphng constant is related to the scat- 
tering amplitude f{E) for Cooper pairs by g~^f = 

This requires us to solve 
the two-body problem for finding the scattering ampli- 
tude in the presence of the ID lattice. In this case the 
expression for f{E) is given by 



f{E) = a / dZcj,*E{Z,0)dr{r'i>{Z,r))r=o 



(5) 



where (j)E{Z,r) = (22)6*'^^'"^ is the incom- 

ing wavefunction for two atoms undergoing Cooper pair- 
ing. The center of mass and relative coordinates are Z — 
{zi +Z2)/2 and r = ri — r2, and E = ?i^q^/m + 2ei(q^) is 
the total energy. The two-particle wavefunction ^'(Z, r) 
obeys the Schrodinger equation 

-A-—— + V{Z,z)+gS{r)-r-EU^O, (6) 
m 4m oZ"^ or J 

where V{Z,z) = Vopt{zi) + Vopt{z2)- The solution of 
Eq.(|SJl can be written as 

^{Z,Y)^4>E{Z,v)+g J dZGEdr'ir'-^ir',Z'))r'=o (7) 

where Ge{j^, Z;0, Z') is the Green function of Eq.© 
with g = 0. The behaviour of the Green function at 
short distances r is governed by the Laplacian term in 
Eq.p yielding G'£;(r, Z; 0, Z') = -5{Z - Z')m/^TTh\ + 
Ke{Z, Z'), where Ke{Z, Z') is a regular function. Then, 
from Eq.((2|) we immediately obtain an equation for the 
function Y{Z) = dr{r'i>(r, Z))r=o appearing in Eq.©: 



Y{Z)^(^E{Z,0)+g J dZ'KE{Z, Z')Y{Z') 



(8) 



Writing the kernel of the integral equation ((HJ in the form 

Ke{Z,Z') = [GB(r,Z;0,Z') - Gi5=o(r,^;0,Z')].=o + 
Ke=q(Z^ Z'), we expand the Green function G^; in eigen- 
statcs of non-interacting atoms: 



(Z') 

E + iO- En^qiz) - en2{q2z) ~ h'^q'i/m 



(9) 



and retain only the contribution of the lowest Bloch band. 
In the tight binding limit, these states can be written in 



terms of Wannier functions as <piq^ (z) ~ J2i 



w{z — 



Id), where w{z) = (l/Tr^/^'cr^/^) exp(— z^/2(t^) is a vari- 
ational Gaussian ansatz. By minimizing the energy of 



non-interacting lattice atoms with respect to cr, one finds 
d/cr = 7rsi/4exp(-l/4Vi) 0. 

We now insert the ansatz Y{Z) — AJ2(U)^{Z — (-d) 
into Eq.Q and take into account that the relation 
^ dZ'KE=o{,Z,Z')Y{Z')dZ' = y(Z)m/47rfi2a gives a 
critical value of the scattering length a = needed to 
form a two-body bound state in the lattice [lO|. Then, 
using the dispersion relation ei{qz) — 2t{l — cos{qzd)) and 
obtaining the kernel Ke{Z, Z') on the basis of Eq.JHJ, we 
find the coefficient A in the expression for Y{Z). Equa- 
tion then leads to the scattering amplitude 



fiE) = 



iC 



1 - a/acr + {a/V2^a)a{E/4t) ' 



(10) 



where C—dl\j2/Ka . The function Q:(a;)=zarccos(l — x) for 
x<2, and Q!(a;) = - ln[a;(l + ^l-2/.x)2/2] +i7r for x>2. 

Equation H10() is one of the key results of this paper. It 
shows that the scattering amplitude undergoes a dimen- 
sional crossover as a function of energy. In the anisotropic 
3D regime {E <^ 8t) we have / = aC/(l — a/ocr + 
iaCvEm* /H) where m* is the effective mass at the bot- 
tom of the band. In the quasi-2D regime {E ^ 8t), the 
tunneling between wells is irrelevant and the two atoms 
are in the ground state of an effective harmonic potential 
of frequency ujq — h/ma^ . The scattering amplitude of 
Ea. (|10|) should then reduce to / = f2Dd where /2D = 
(a/\/2^CTl/ (1 + {a/y/2^a){hi[\h^o/E] + iir)) and A = 
0.915/7r jsj. This provides us with the asymptotic be- 
haviour of the critical value of the scattering length 

a^r = -V2TT(j\\-Y^{\huoal2t), (11) 

which agrees with numerics of Ref. 10] already for s ^ 5. 

We see that the ID lattice affects f{E) and the cou- 
pling constant 5e// in a non-trivial way. For e^? < 4t, the 
coupling is density independent and from Eq.(|3Jl we get 



^0 27 _p 

J-c — — 6 EE exp 

TT 



2qzFC 



1 



1 



(12) 



where the function F is given by Eq.Q), and qzF — 
arccos(l — ep/^tj/d is the Fermi wavevector along the 
z-axis. Equation H12I) is valid provided ^ ep, which 
implies \a\ < \acr\- In the low density limit ep <C 4t, 
reduces to the mean field transition temperature [ll| 
for a homogeneous gas of atoms with an anisotropic 
quadratic dispersion and a renormalized 3D inverse scat- 
tering length a 
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C^^{\a\~^ — |acr|^^)- The presence 
of the lattice causes an effective shift of the resonance 
from 1/a = to 1/a = l/a^r < 0, which in turn gives 
rise to a sharp increase in at a fixed value of the 3D 
scattering length. 

For ep > 4:t, the coupling constant ffe// becomes den- 
sity dependent. Equation igj yields F{fi) = 21n[2/(l + 
y/1 ~ 4t/^)], and from Eq.© we find 



(13) 
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Note that the exponent in the rhs of Eq. H13() does not de- 
pend on the Fermi energy, the density of states being con- 
stant for ep > At. For ep ^ At, the ratio At/ep plays 
the role of a small parameter ensuring the inequality 

/ep <^ 1 also for values of the scattering length larger 
than Ocr, but still \a\ ^ cr. In this regime, the system 
behaves as a stack of quasi-2D superfluid gases weakly 
coupled by Josephson junctions. The transition temper- 
ature in each disc takes the form T° = "fy/2epEb/'!T |12l |. 
where in our case Et is the binding energy of the two- 
body bound state in the ID optical lattice. In the absence 
of coupling between the discs {t = 0), the gas in each disc 
is two-dimensional and the transition is therefore of the 
Kosterlitz-Thouless type. However, for 2D BCS super- 
fluids a standard calculation of the Kosterlitz-Thouless 
transition temperature gives a value that is lower than 
Tc by an amount ~ T'^jep <C 0| and thus lies inside 
a narrow region of critical fuctuations in the neighbor- 
hood of Tc- Therefore the mean field approach leads to 
a correct result for the transition temperature. 

We next proceed to evaluate Gorkov's correction to 
the transition temperature due to the polarization of the 
medium Following Ref.|3li we introduce the static 
Lindhard function 



i(p) 



/(fq) - /(Cq+p) 



(27r)3 



(14) 



where j(x) = Q{—x) is the Fermi distribution at T = 
0. The induced interaction between two states q and 
q' on the Fermi surface is given by Vindi^') = 9CL{p), 
with p = q -I- q' and C defined as above. Since L{p = 
0) = i^ifJ,), we write L{p) = v{^)B{p), where _B is a 
dimensionless positive function sensitive to the geometry 
of the Fermi surface. The critical temperature is then 
given by = T^e'^^^'"', where 

^""^"^ - — Jsium^^' — ■ ^ 

The integration in (|15|) is done numerically and the cor- 
responding Gorkov correction Tc/T^ is shown in Fig^ 
For ep <C At, the system has an anisotropic quadratic 
dispersion and we recover the result for the homogeneous 
case {B)ps = (1 -h 21n2)/3, yielding = T°/{Ae)'^/^ = 
0.45r". In the hmit ep » At, the band dispersion e(fc2) 
can be neglected in Ea. lfTHl and we find Tc/T^ — e^^ , in 
agreement with Ref.[l|. The cusp at ej? = 4t is expected 
as this is the point of the Van Hove singularity 0| where 
the derivative of the density of states, dv/de, diverges. 

So far we have discussed the BCS superfluid transition 
in a ID optical lattice. In the second part of the Letter 
we assume that the superfluid gas is at zero temperature 
and it is confined in the x, y directions by a trapping 
potential. Then, as the tunneling rate between neighbor- 
ing discs is tuned below a critical value tc, the system 
undergoes the superfluid-Mott insulator quantum tran- 
sition. For a large number of atoms per well {N 3> 1), 
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FIG. 1: Gorkov's correction versus ep/At. The limiting value 
Tc/T° = e"^ at ep/At > 1 is shown by a dotted line. 



the critical hopping rate can be evaluated within the hy- 
drodynamic approach . Neglecting the coupling with 
radial degrees of freedom and the particle loss due to in- 
elastic processes, the proper dynamical variables are the 
particle number fluctuation N'^^ and the phase of the 
order parameter in each disc. The hydrodynamic equa- 
tions are equivalent to the classical equations of motion 
of the ID phase Hamiltonian 

Hp = J2iEc/2)Nf ~ Ej cos($,+i - $£), (16) 



where Ec — and Ej = t^N / ^ are the charging and 

the Josephson energies, respectively, and HN'^ and $f are 
considered as conjugated variables. 

Quantization of the classical Hamiltonian H16|) is 
achieved by replacing these variables with operators hN' 
and $ satisfying the commutation relation [hN' , = ih. 
The quantized Hamiltonian is known to exhibit a phase 
transition at the critical value Ec = TjEj, with 77 ~ 0.81 
The superfluid phase occurs for Ec < tjEj and is 
characterized by an algebraic decay of the phase correla- 
tion function (cos($£ — $fe)) at large distances \i—k\ ^ 1. 
The decay becomes exponential for Ec > rjEj where one 
enters the Mott regime, characterized by large phase fluc- 
tuations which suppress interwell tunneling. The ground 
state is an insulator with a fixed number of atoms per disc 
and a finite gap in the excitation spectrum. By compar- 
ing the values of the charging and the Josephson energies, 
we find that for BCS superfluids 



t_c 



1 n 



(17) 



This result differs from the corresponding value for Bose 
condensates in the same geometry, t^//^'' ~ where 
IJp = ng2D is the chemical potential and g2D is the 2D 
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coupling constant. This is because the Josephson energy 
in the Hamiltonian (|16|l for the bosonic case is = tN. 

Equation (|17|l has been derived under the assumption 
that the effective tunneling rate i^c = tc/^l^ is large com- 
pared to the loss rate j>. The most severe losses come 
from three-body recombination. For an array of Bose- 
condensed atomic gases in the same geometry the corre- 
sponding loss rate is always large compared to the crit- 
ical tunneling rate 0, unless the number of atoms per 
disc is very small as in the experiment of Ref. 0. For 
Fermi superfluids the situation is completely different be- 
cause the inelastic processes are strongly inhibited by 
quantum statistics. In the quasi2D geometry, in anal- 
ogy with the 3D case 0|, for the 3-body loss rate one 
can write /> = —h/n = Lri^{kFRe)^ , where the small 
factor (kpRe)^ comes from the Pauli principle, Re is a 
characteristic radius of the interatomic potential, and L 
is the quasi2D recombination coefficient. This coefficient 
is related to the corresponding quantity of a 3D gas as 
L ~ L^jj/a^. Within an order of magnitude, L^jj co- 
incides with the recombination rate constant for bosonic 
isotopes of the same atom, which ranges from 10^^^ to 
10~'^''cm^/s for alkali atoms. The ratio of the loss to the 
effective tunneling rate is then given by 



2mL, 



3D 



(18) 



and should be much smaller than unity for consistency. 

Note that in the case of bosons the effective tunnel- 
ing rate is — t^/h ~ 'ng2D/N'^- For a given iV, 
it is smaller than the one for fermions by a factor of 
{fTig2D / 1^^^) , which is a small parameter of the the- 
ory for the 2D weakly interacting gas. The bosonic 3- 
body loss rate is ^ (i3£)/cr^)n^ and it exceeds the 
fermionic loss rate by a factor of (kpRe)~'^. Thus, we 
have iy/vc ~ {'rag2DlT^^^){kFRe)'^v^ /v^- For realistic 
densities the fermionic ratio v/vc is smaller than the 
bosonic ratio lv\ by 4 or 5 orders of magnitude. This 
is a consequence of Fermi statistics and it is crucial for 
the observation of the Mott transition. 

For example, considering N — 10^ fermionic potas- 
sium atoms (i?e K, 5 nm) in each disc, with density 
n = lO^cm^^ corresponding to a Fermi energy (chem- 
ical potential) tp — [i = 380 nK k^, from Ea. H17|) we 
find ic/ft ~ 70 s^^ corresponding to s 25 for a lat- 
tice period d — 400 nm, and cr w 60 nm. This leads to 
Vc ~ 0.1 s~^ and, assuming < 10~^* cm^/s, from 

Ea. p8|l we obtain v jvc < 0.1. Hence, owing to quantum 
statistics, in Fermi superfluids one can easily have v <^Vc 
and achieve the Mott insulator transition. It is important 
to emphasize that, in the given example, the suppression 
of recombination processes by a factor of {kpRe)"^ ~ lO"'^ 
originating from the Pauli principle is crucial to keep the 



ratio (EHJ small even for N 10^. 

In conclusion, we have found the superfluid transition 
temperature for a two-component Fermi gas in a ID op- 
tical lattice and revealed that the effective coupling con- 
stant depends in a non-trivial way on both the atom den- 
sity and the parameters of the optical field. For an array 
of finite size discs with a large number of atoms per disc, 
we have shown that the critical effective tunneling rate 
for the Mott insulator quantum transition can be larger 
than the rate of particles losses. Thus, the Mott phase 
can be observed for Fermi superfluids in this geometry. 
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